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Abstract. In simulations of some economic gas-like models, the asymptotic regime shows an exponential 

of this statistical equilibrium for this type of gas-like models is explained in a rigorous analytical way. 

^^—^ 1 Introduction duce the two different types of statistical behavior before 

I— I mentioned. The exponential distribution can be obtained 

Exponential (Gibbs) distribution is ubiquitous in many , . . , , . , . 

by supposmg a gas oi agents that exchange money m bi- 

CO natural phenomena. It is found in many equilibrium situ- .... . „ .... . . , 

nary collisions, or m hrst-neighbor interactions, with ran- 

ations ranging from physics to economy. In an economic . , ... ... ,. . , , . 

l/^ dom, deterministic or chaotic ingredients m the selection 

context, the society of the Western (capitaUst) economies ^ , . , , ■ ™ 

^— 5 of agents or m the exchange law |4l5l6j . ihe power law 

behavior is obtained by introducing some kind of inho- 



o 

can be divided into two distinct groups according to the in- 



comes of people [T]. The 95% of the population, the middle 



mogeneity in the gas through the saving propensity of the 



and lower economic classes of society, shows an exponen- _ , , , , , . r ^ 

5-H agents 7 or through the breaking of the pairing symmetry 

^ . 

tial wealth distribution. The rest of the population, the . , , , r^i 

m the exchange rules [6j . 

5% of individuals, arranges their privileged incomes in a 

power law distribution. In this work, we take a continuous version of one of 

Gas-like models interpret economic transactions be- such models [5] and we show that the exponential distri- 

tween agents trading with money similarly to collisions bution is the only fixed point of the discrete time evolution 



in a gas where particles exchange energy |2|3| . Different operator associated to this particular model. This rigor- 
gas-like models have successfully been proposed to repro- ous analytical result helps us to understand the computa- 



2 Title: Exponential wealth distribution in 

tional statistical behavior of such a model which reaches 
its asymptotic equilibrium on the exponential distribution 
independently of the out of equilibrium initial state of the 
system. 

2 The continuous gas-like model 

We take the simplest gas-like model [1] in which an en- 
semble of economic agents trade with money in a ran- 
dom manner. An initial amount of money is given to each 
agent, for instance the same to each one. Transactions be- 
tween pairs of agents are allowed in such a way that both 
agents in each interacting pair are randomly chosen. Also 
they share their money in a random way. When the gas 
evolves with this type of interactions, the exponential dis- 
tribution appears as the asymptotic wealth distribution. 
In this model, the local interactions conserve the money, 
then the global dynamics is also conservative and the total 
amount of money is constant in time. 

The discrete version of this model is as follows 3!- 
The trading rules for each interacting pair (m^, rrij) of the 
ensemble of N economic agents can be written as 

m'i — e {rrii -f nij), 

m'j = (1 - (){mi +mj), (1) 
i,j = 1...N, 

where e is a random number in the interval (0,1). The 
agents (i, j) are randomly chosen. Their initial money (m^, r 
at time t, is transformed after the interaction in (m,^, m^ ) 
at time t+1. The asymptotic distribution p^(to), obtained 
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by numerical simulations, is the exponential (Boltzmann- 
Gibbs) distribution, 

P/(m) = /3 exp(— /3 m), with P = l/<m>gas, (2) 

where p/(m) dm denotes the PDF {probability density func- 
tion), i.e. the probability of finding an agent with money 
(or energy in a gas system) between m and m + dm. Evi- 
dently, this PDF is normalized, ||p/|| — Pf{m)dm — 1. 
The mean value of the wealth, < m >gas, can be easily 
calculated directly from the gas by < to >gas= X^i Tfii/N . 

The continuous version of this model consists in think- 
ing that an initial wealth distribution po^m), with a mean 
wealth value < po >, evolves from time n to time n + 1 
under the action of an operator T to asymptotically reach 
the equilibrium distribution pf{m), i.e. 

lim T" (po (m) ) ^ (to) . (3) 

In this particular case, Pf{m) is the exponential distribu- 
tion with the same average value < pf >=< po >, due 
to the conservation of the total money. This is a macro- 
scopic interpretation of Eqs. ([T]) in the sense that, under 
this point of view, each iteration of the operator T means 
that many interactions, of order N/2, have taken place 
between different pairs of agents. As n indicates the time 
evolution of T, we can roughly assume that t ^ N n/2, 
where t follows the microscopic evolution of the individ- 
ual transactions (or collisions) between the agents (this 
alternative microscopic interpretation can be seen in [3]). 
j), Now, we recall how to derive T (see Ref. [8]). Sup- 
pose that p„ is the wealth distribution in the ensemble at 
time n. The probability that two agents with money (m, v) 
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TO OBTAIN X 




Prob{{u,v) X) 



p(u)(Ju ' p(v)clv 



Fig. 1. Derivation of the nonlinear evolution operator T 
for the agent-based economic model given by Eqs. ([T]). See 
explanation in the text. 

interact will be Pn{u) *pn{v). As the trading is totally ran- 
dom, their exchange can give rise with equal probability 
to any value x comprised in the interval (0, u + v). Then, 
the probability to obtain a particular x (with x < u + v) 
for the interacting pair (u, w) will he Pn{u) *pn{v)/(u + v). 
Finally, we can obtain the probability to have money x at 
time n + 1. It will be the sum of the probabilities for all 
the pairs of agents (u, v) able to generate the quantity x, 
that is, all the pairs verifying u + v > x. T has then the 
form of a nonlinear integral operator (see Fig. [T]) , 

Pn{u)pn{v) 



Pn+l{x) = Tpn{x) 



U + V 



dudv . (4) 



If we suppose T acting in the PDFs space, we prove in 
the next section that T conserves the mean wealth of the 
system, < Tp >=< p >. It also conserves the norm (|| • ||), 
i.e. T maintains the total number of agents of the system, 
llTpll — \\p\\ = 1, that by extension implies the conser- 
vation of the total richness of the system. We also show 
that the exponential distribution Pf{x) with the right av- 



erage value is the only steady state of T, i.e. Tpf = Pf, 
and that order two periodic points do not exist. It also 
seems feasible that other high period orbits do not exist. 
In consequence, it can be argued that the relation Q is 
true. See for instance simulations done in Fig. [2] The same 
result is obtained for any other initial distribution [TU] . 




Fig. 2. An example of the action of T in the space of 
PDFs. In this case, po{x) is chosen as a normalized trian- 
gular distribution. Observe that after two iterations of T 
the system has practically reached the asymptotic state, 
i.e. the exponential distribution Pf{x). (Figure avanced 
from M)- 



3 Properties of operator T 
3.1 Definitions 

Definition 1 We introduce the space of positive func- 
tions (wealth distributions) in the interval [0,cxd), 

L+[0,oo) = {y: [O,oo)^i?+U{0}, ||y|| <oo}, 

with norm 

112/11 = / y{x)dx. 
Jo 
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Definition 2 We define the mean richness < x >y asso- 
ciated to a wealth distribution y G ij^[0,oo) as the mean 
value of X for the distribution y. In the rest of the paper, 
we will represent it by < y >. Then, 

1*00 

< y >=< X >y=\\xy{x)\\^ I xy{x)dx. 

Jo 

Definition 3 For x > and y E Lf[Q,oo) the action of 
operator T on y is defined by 

y{u)y{v) 



Tivix)) 



dudv- 



IS(X) U + V 

with S{x) the region of the plane representing the pairs of 
agents {u, v) which can generate a richness x after their 
trading, i.e. 

S{x) = {(u, v), u,v > 0, u + V > x}. 

After some manipulations, other representations of the 
integral operator T are possible: 



{Ty){x) = X 
{Ty){x) = X 
{Ty){x) = 



1 nOO poo POO 

dt ds+ dt ds 
Ji-t Ji Jo 

dt / ds- dt ds 
"'0 "'0 Jo 



y{xt)y{xs) 
t + s ' 

y{xt)y{xs) 
t + s ' 



dr / du y(ur)y{{l—u)r). 

X Jo 



3.2 Theorems 

Theorem 1 For any y <E L^[0,oo) we have that \ \Ty\ \ = 
In particular, for y being a PDF, i.e. if \\y\\ — 1, 
then \\Ty\\ — 1. (It means that the number of agents in 
the economic system is conserved in time). 



Proof Take y G L{[0,oo). Then (see Fig.|3 

yiu)y{v) 



\\Ty\\ = / dx 
Jo 



dudv- 



S(x) 



U + V 



^du Hdv rdxy^ 

Jo Jo u + v 



00 



CO 



y{u)du / y(v)dv — \\y\ 
^0 



Corollary 1 It is clear that {Ty){x) > for x > and 
Ty G i]^[0, 00). Also, it is straightforward to see that Ty is 
a continuous function, i.e. Ty G C[0, 00). Therefore, from 
the above Theorem we have that Ty G L'^lQ, 00) p| C[0, 00), 
that is T : L'l[Q,oo) — > ij^ [0, 00) f] C[0, cx)). 

Corollary 2 Consider the subset of PDFs in ij^[0, 00), 
i.e. the unit sphere B := {y G L]^[0, 00), \\y\\ = 1}. Ob- 
serve that if y (z B then Ty G B. Therefore, we can 
define the restriction Tb of T to the subset B , that is 
Tb-.B^B. 



(0,x,x) 




Fig. 3. Integration region in the triple integral considered 
in the proof of Theorem [T] 

Proposition 1 The operator T is Lipchitz continuous in 
B with Lipchitz constant < 2. 
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Proof. Take y,w G B. Then 



\\Ty-Tw\ 



Jo Jo u + v 

oo poo 



du / dv\y{u)y{v) — y{u)w{v) + 
Jo 

+y(u)w{v) — w{u)w{v)\ < 

POO poo 

< / y{u)du / dv\y{v) — w{v)\ + 
Jo Jo 

/•oo /'OO 

+ / w{v)dv / du\y{u) — w{u)\ = 2\\y — w\\. ^ 
Jo Jo 

Theorem 2 The mean value < y > of a PDF y is con- 
served, that is < Ty >—< y > for any y Cz B. (It means 
that the mean wealth, and by extension the total richness, 
of the economic system are preserved in time). 



Proof. 

<Ty>= \\x{Ty)ix)\\ = I xdx 







dtds^^^ 

S(x) t + S 



^dt Hds rxdx^^ 

Jo Jo t + S 

poo poo 



dt / ds{t + s)y{t)y{s) 



'0 Jo 

^ poo 



ty(t)dt + - sy{s)ds = 
Jo ^ Jo 

POO 

= / xy{x)dx = \\xy{x)\\ ^< y > . ^ 
Jo 

Observation: Because \\Ty\\ — ||y|p in L^[0, oo), this 
operator divides the space [0,oo) in three regions: the 
interior of the unit ball Bin — {y i^[0, oo), ||?;|| < 1}, 
the unit sphere B defined in Corollary 2, and the exte- 
rior of the unit ball B^^ = {y G L^[0,oo), \\y\\ > 1}. 
Accordingly, the recursion 

yn{x) = T'^yoix), with ya^x) G Lf[0,oo), 

with the norm 1 1 • 1 1 has three different behaviors depending 
on yo: 



- ff I/O e Bin, then lini„^oo ||yn(a;)|| =0. 

^ If 2/0 G B, then ||?/n(2^)|| = 1 (the succession yn remains 
in B y n). 

- If yo e B^^, then lim„_^oo l|yn(a;)|| = oo. 

Then, it is obvious that y = is a fixed point of T. 
The remaining fixed points of T in -/j^[0, oo), if any, must 
be located in B. In fact, we have the following Theorem. 

Theorem 3 Apart from y — 0, the one-parameter family 
of functions ya{x) — ae^"^ , a > 0, are the unique fixed 
points of T in the space Lj^[0, oo). 

Proof. They are fixed points. Using the representation ([s]) 
we find (see Figure 111) that 



2 poo pr -ar 



ale "^dr = ae 



To see the uniqueness, consider the Fourier transform 
of y G L^[0, oo): 



noo 

y{p) = / e'P^y{x)dx, 
Jo 



and the Fourier transform of Ty G if [0, oo): 

{T^){p)^ r e^P^dx f f ^y{u)y{v) 
Jo J Js(x) u + v 



CO POO pu+v ^ipx 

duy{u) / dvy{v) / dx 

Jo Jo u + v 



1 



duy{u) I dvy{v)- 
Jo Jo ^ + 

If {Ty){x) = y{x) then {Ty){p) = y{p) and 



poG POO ^ip(u+v) ^ 

ipy{p) = / duy{u) / dvy{v) . 

Jo Jo u + v 

Taking the derivative with respect to p in both sides of this 

equality we find that the Fourier transform y{p) satisfies 
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the differential equation: 

y{p) +py'{p) = y'^ip), 

whose general solution is 



yip) 



c€C. 



1 ~ cp 

When we are considering only positive functions y{x), 
their Fourier transform must satisfy y{p)* — y{^p)- This 
is only possible if c = — c*, that is, the constant c must be 
purely imaginary: c ^ i/a, a Cz R and then 



y{p) 



a e R. 



1 - {i/a)p' 

From the inverse theorem for the Fourier transform [[II , 
p. 61, Theorem 8.3] we have that any solution y{x) of the 
equation Ty = y in i]^[0, oo) is of the form: 



1 



oo ^-ipx-e\p\ 



y{x)^ — hm / dp. 

2ti e->oy_^ 1 - {i/a)p 

From [[H], p. 197, Lemma 1] we have that 



1 

2^ 



-dp. 



1 - {i/a)p 

Because x > 0, the integrand is exponentially small when 
the imaginary part of p is large and negative. Then, we 
can close the integration line (— oo, oo) with a semi-infinite 
circle on the lower complex p-plane: 



y{x) 



-dp, 



27r Jc 1 - 

where C is the union of the real line and that semicircle 
traveled in the negative sense (clockwise). The pole p — 
—ia of the integrand is outside of C if a < and inside 
C if a > 0. Applying the Cauchy Residua Theorem to the 
above integral we obtain that 

if a < 0, 

ae""^ if a > 0. 



yix) 



Then, y = and ya{x) — ae a > 0, are the unique 
fixed points of T in the space ^^^[0, oo). ^ 



VA 




U 



Fig. 4. Integration region in the double integral consid- 
ered in the proof of Theorem |3j 



Corollary 3 The family of functions ya{x) = ae^°"^ , 
a > 0, are the unique fixed points of Tg (the restriction 
of T in the unit hall B). 

As we will show with some examples in the next sec- 
tion, the numerical results indicate that, apart from the 
family of exponential distributions, there is no other peri- 
odic or aperiodic asymptotic sets for the operator T acting 
on the PDFs ball B. Thus, the whole unit ball B seems 
to contract to the exponential steady-state regime. In the 
next Theorem, we make another step in this direction by 
showing that the operator T has no cycles of order two. 

Theorem 4 The operator T has not cycles of order two 
in L]^[0, oo), that is, T^y = y Ty = y. 



Title: Exponential wealth distribution in a random market. A rigorous explanation 



Proof. A similar computation to the one introduced in the 
proof of Theorem |3] shows that 

{Ty){p) ■ ' - ' ' - 



(_l)™(T"y)M(0) 



m—l 



-^(r"-i2;)W(o)(r"-iy) 

m ^ — ^ 



{m-k-l) 



(0). 



OO nOO 

du{Ty){u) / dv{Ty){v)— — ; — — 
W(u + v) 



fc=0 



v{p)- 



Proof. For fixed n, let us show that {T"y) is m times con- 
tinuously differentiable and e L+[0,oo) 
by induction over m. It is obvious for m = 0: {T^y) is 
continuous and (T^y) G -^^[0, oo) for any n — 1, 2, 3, . . . 



Taking the derivative with respect to p in both sides of 

both equations, we obtain the foUowinff system of differ- c ^ ■ i- r ^ ■ i^ jw 

^ ' " buppose that it is true tor a certain m — 1 < n and let 



us 



ential equations for y{p) and [Ty){p): 
{[p{Ty){p)]' ^ tip). 

\ [py{p)]'^ m\p)- 

Taking the cross product of these equations and multiply- 
ing by p^ we get the equation: 

P^t{p)[pv{p)]' - p\Tyf{p)[p{T^){p)]' , 
and integrating: 



iTy)Hp)^yHp) + 



c e 



But both, y{p), {Ty){p) € C(— oo, oo) and then c = 0. This 
means {Ty){p) = y{p) and {Ty){x) = y{x) (a.e.). ^ 

To finish this section, another interesting property of 
operator T concerning to the value of the norm of the 
derivatives of r"y is disclosed in the following Proposition. 
A consequence of this property is also given. 

Proposition 2 For any y e Lj^[0, oo) and m,n G N, 

with m < n, we have the following properties: 

{T^y) is m times differentiable, 

(-l)"(T"y)('")eL+[0,(X3), 
ll(T"2/)(")|| = |(T"2/)("-i)(0)|, 



show that it is also true for to. From the last representa- 
tion of the integral operator T given after Definition [3| we 
have that 



{r^yYix) = - [\T-~^y){ux){r^-\j){{l - u)x)d 
Jo 



u < 0, 



and then 



(-i)™(r"y)(™)(x) 



TO — 1^ ' 



E 

fc=0 



\-l)'^{T--^y)i'')(ux) 



m—k — 1 /rp'fi—l \{7n—k—l) 



i{l-u)x) 



du. 



From the induction hypothesis, all the integrands in the 
above sum of integrals are continuous and then (T"?;)*^™' 
is continuous and integrable. Moreover, they are positive 
and then (-1)"(T"?;)(") (x) > 0. The last two assertions 
of the Proposition are trivial. ^ 

As a consequence of the former Proposition, we have: 

Proposition 3 A fixed point y of T , if any, must be a 
completely monotonic function: y G C°°[0, oo) and (— l)™?/'™^ (x) > 
Vx £ [0, oo). Then, as it can be easily tested, the expo- 
nential distribution ya{x) = ae~°'^ , a> 0, is one of such 
type of functions. 
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Proof. If Ty ~ y for a certain y £ L^[0, oo), then r"y = y We numerically check that 
for any n £ N. For arbitrary m £ N , take n > m. Our 
goal is proved just seeing from the former Proposition that 
(-l)'"y('")(.T) = (-l)™(T"j/)(™)(j;) > 0. J|k 



\\Tya,n -W\\< WVam " w\ 

for different values of a and n. 



4 Time evolution of T: Examples 



Example 2 For a, /3 > 0, consider the new family 



yc.A^) = l[ae-"^ + l3e-^-] £ B. 



The mean value is: 

a/3 



As discussed in the former section, the iteration of T for 
any PDF y in the unit ball B with a finite < y > seems to 

converge toward the only fixed point of the system, i.e. the < Va p 7^^ 

exponential distribution presenting the same mean value 
of the initial condition y. In this respect, we make the 
following conjecture: 

Conjecture For any y £ B with a finite < y >, {Tya,i3){x) 



2(a + /3)' 

The limit distribution should be w{x) — 76^^^. The first 
iteration of T gives 



with 



— =< y >= / xy{x)dx. 



Again, we numerically check that 

\\Tya,l3 -W\\< \\ya,li - W\\. 

From simulations (see Fig. [2]), it can be guessed that 

the convergence suggested in this conjecture is fast. Next, ^^"^ several values of a and /?. 

we compute and verify the contraction of the dynamics e^^^^i^ 5 For « > and n G N, consider the e-family of 

toward the exponential distribution under the action of T functions 
for some particular PDFs. 



y.a,„(x) = (1 - e)ae-"^ + e^x"e-«^ 



Example 1 For a > and n G N, consider the family 



Q/"+i _ with < e < 1. It is straightforward to see that y^ a n ^ B. 



The mean value is: 



The mean value is: 

_i _ 1 + en 



1 , , , < y^,oc,n >= 7 

7" - ('^ +!)/«■ " 

The limit distribution of T"y„,„ when n ^ 00 should be ^he limit function should be w{x) = 76"^^ The first 

w{x) = 76-'^^. The first iteration is: iteration is: 

(rya,„)(a;) = ^2n+i^|^^ ^ 3/2) ^^^^ + ^' T^^e.^.^l^) = a |l + e 
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